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Abstract 

This paper presents the design of the X-33 Flush 
Airdata Sensing (FADS) system. The X-33 FADS uses a 
matrix of pressure orifices on the vehicle nose to 
estimate airdata parameters. The system is designed with 
dual-redundant measurement hardware, which produces 
two independent measurement paths. Airdata parameters 
that correspond to the measurement path with the 
minimum fit error are selected as the output values. This 
method enables a single sensor failure to occur with 
minimal degrading of the system performance. The 
paper shows the X-33 FADS architecture, derives the 
estimating algorithms, and demonstrates a mathematical 
analysis of the FADS system stability. Preliminary 
aerodynamic calibrations are also presented here. The 
calibration parameters, the position error coefficient (e), 
and flow correction terms for the angle of attack (5a), 
and angle of sideslip (5(3) are derived from wind tunnel 
data. Statistical accuracy of the calibration is evaluated 
by comparing the wind tunnel reference conditions to the 
airdata parameters estimated. This comparison is 
accomplished by applying the calibrated FADS 
algorithm to the sensed wind tunnel pressures. When the 
resulting accuracy estimates are compared to accuracy 
requirements for the X-33 airdata, the FADS system 
meets these requirements. 

Nomenclature 

a flush airdata sensing geometry 

coefficient 

A angle of attack triples algorithm 

coefficient 
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A/D analog-to-digital conversion 

A 1 first angle of sideslip triples algorithm 

coefficient 

^(Mo©) zero’th order angle of attack 

calibration coefficient 

A [(M^) first order angle of attack calibration 

coefficient 

^ 2 (^ 00 ) second order angle of attack 

calibration coefficient 

A 3 (M 00 ) third order angle of attack calibration 

coefficient 

b flush airdata sensing geometry 

dummy variable 

B second angle of sideslip triples 

algorithm coefficient 

B' first angle of sideslip triples algorithm 

dummy variable 

) zero’th order angle of sideslip 

calibration coefficient 

^{(M^ ) first order angle of sideslip calibration 

coefficient 

# 2(^00 ) second order angle of sideslip 

calibration coefficient 

^(M^ ) third order angle of sideslip calibration 

coefficient 

CTE coefficient of thermal expansion 

C second angle of sideslip triples 

algorithm dummy variable 

DFRC Dryden Flight Research Center 

ESP electronically scanned pressure 

module 

/[M^ ] subsonic subsonic Mach number function 
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/[M oolsupersonic 

supersonic Mach number function 

5P 

angle of sideslip correction term, deg 

FADS 

flush airdata sensing 

r 

pressure difference for triples 

Fail Op 

fail-operational (airplane operates as if 


algorithm, lb/ft 2 

there were no failures) 

Y 

ratio of specific heats 

M 

estimation algorithm geometry matrix 

E 

position error calibration parameter 

MAAF 

Michael Army Air Field, Tooele, Utah 

e M 

variation of E with Mach number at 

MAFB 

Malmstrom Air Force Base, Great 
Falls, Montana 

e al 

a f , = 0° 

first order coefficient, fit of E with 

MSE 

mean square error 


angle of attack 

second order coefficient, fit of £ with 


free stream Mach number 

a 2 

angle of attack 

NASA 

National Aeronautics and Space 

e pi 

e [12 

first order coefficient, fit of e with 

P 

Administration 
port pressure, lb/fr 

angle of sideslip 

second order coefficient, fit of £ with 

Poo 

2 

free stream static pressure, lb/ft 

angle of sideslip 

Q 

weights on measured pressures 

X 

cone angle of FADS port, deg 

Qc 

impact pressure, lb/fr 

<t> 

clock angle of FADS port, deg 

Q 

pressure weighting matrix 

e 

local flow incidence angle, deg 

RCC 

reinforced carbon-carbon 

a 

standard deviation 

RLV 

Reusable Launch Vehicle 

y 

dummy variable used in angle of 

RMS 

root mean square 


sideslip root analysis 

TPS 

thermal protection system 

x" 

chi-squared distribution 

2 . 2 

UPWT 

Unitary Plan Wind Tunnel, Langley 
Research Center, Hampton, Virginia 

Q FADS geometry term, cos U + esin u 

Superscripts, Subscripts, and Mathematical Operators 

X 

state vector 

i 

port index 

X 2 

2 2 

mean square fit error, (lb/ft ) 

j 

port index 


a 

a e 
® true 

P 

P, 

P true 

Pi 

P2 

5a 


angle of attack, generic, deg 

local angle of attack sensed by FADS, 
deg 

free stream angle of attack, deg 

angle of sideslip, generic, deg 

local angle of sideslip sensed by 
FADS, deg 

free stream angle of sideslip, deg 

quadratic solution 1 for angle of 
sideslip, deg 

quadratic solution 2 for angle of 
sideslip, deg 

angle of attack correction term, deg 


(]) 

k 

A 

Det[.] 

I .]' 1 

if 

V[.] t 

3M 


iteration index 

port index 

FADS estimate 

matrix determinant 

matrix inverse 

matrix transpose 

gradient with respect to vector x 

partial derivative with respect to x 

Introduction 


The primary goal of the Single-Stage-to-Orbit 
Technology program is to radically reduce the cost of 
access to space, and the X-33 advanced technology 
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demonstrator is the centerpiece of this effort. The X-33 
design is a 53-percent scale model of the Lockheed 
Martin VentureStar™* Reusable Launch Vehicle (RLV) 
and is based on a lifting-body design which features two 
linear aerospike rocket engines. The autonomous-flight 
X-33 will launch from Edwards Air Force Base at 
Edwards, California and land at one of two sites: 
Michael Army Air Field (MAAF) in Tooele, Utah, or 
Malmstrom Air Force Base (MAFB) in Great Falls, 
Montana. The vehicle is designed to achieve a peak 
altitude near 300,000 ft and speeds of greater than Mach 
12. After atmospheric re-entry, the X-33 returns to Earth 
with an unpowered horizontal landing. A comparison of 
the X-33 and VentureStar™ vehicles is presented in 
figure 1. 

Because the X-33 is to perform an unpowered landing, 
knowledge of the dynamic pressure, angle of attack, and 
surface winds is critical, so that the terminal area energy 
management (TAEM) can ensure that the target runway 
be reached under a wide variety of atmospheric 
conditions. Direct feedback of angle of attack and angle 
of sideslip are also required for gust load alleviation on 
the vehicle airframe during the ascent of the flight. Thus, 
it was determined early in the X-33 program that the full 
airdata state, including Mach number, angle of attack, 
angle of sideslip, dynamic pressure, airspeed and 
altitude, would be a flight-critical requirement for both 
the RLV and the X-33. 

To achieve the TAEM and ascent airdata requirements 
the airdata system must be operational between 
Mach 0.20 and Mach 4.0 and meet the following 1-a 
accuracy requirements. 

1. Mach Number: ±5.0 percent accuracy for 

2.5 < <4.0, ±2.50 percent accuracy for 

0.6 < M^ <2.5, ±0.015 absolute error between 
M^ = 0.20 and M 00 = 0.60, 

2. Angle of attack : ±1.5° absolute accuracy for first 
three flights, ±0.50° thereafter, 

3. Angle of sideslip: ±0.5° absolute accuracy for all 
flights. 

4. Geopotential altitude: ±200 ft absolute accuracy 
for 0.2 < M^ < 4.0. 

5. Dynamic Pressure: ±15 lb/ft 2 for 0.2 < M^ <4.0. 

To determine the best means of meeting the airdata 
requirements, NASA Dryden Flight Research Center 
(DFRC) performed a feasibility study to compare the 
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performance and cost of a Flush Airdata Sensing (FADS) 
system to a set of deployable probes similar to that 
system installed on the space shuttle. This study 
concluded that a FADS system was more economical 
by a factor of approximately two. Two issues 
made the probe-based system prohibitively expensive: 
1) integration onto the X-33 airframe, and 2) system 
calibration. The FADS system requires no deployment 
mechanisms and can be integrated directly onto the 
vehicle nosecap with no movable parts. Because the 
FADS system does not probe the flowfield, but instead 
uses the natural contours of the forebody, the flow field is 
much cleaner and is easier to calibrate. An additional 
advantage of the FADS system is that it offers the 
potential to sense airdata on ascent, an option not 
available to the probe-based system. Based on the results 
of this study, Lockheed Martin Skunkworks, Palmdale, 
California, selected the FADS system in favor of the 
deployable probes. 

The DFRC FADS design builds on work which 
originated in the early 1960's with the X-15 program, 1 
continued at NASA Langley, 2 ’ 3 and Dryden Flight 
Research Centers 4, 5 in the 1970’s and 1980’s, and 
recently concluded flight testing of an onboard real-time 
system in the early 1990’s. 6, 7 The FADS concept, in 
which airdata are inferred from nonintrusive surface 
pressure measurements, does not require probing of the 
local flow field to compute airdata parameters. This 
innovation allows the extreme hypersonic heating 
caused by the small radius of a flow-sensing probe to be 
avoided, which extends the useful range of the airdata 
measurement system to the hypersonic flow regime. This 
paper describes the X-33 FADS system design, develops 
the aerodynamic model which relates the airdata 
parameters to the measured pressures, and derives the 
computational algorithms used to compute the airdata 
from the pressure measurements. Wind tunnel 
calibration and validation of the FADS system is also 
presented. 

Flush Airdata Sensing Aerodynamic Model 

The fundamental concept of the FADS system is that 
airdata parameters can be estimated from flush surface 
pressure measurements. To perform this estimation, the 
airdata states must be related to the surface pressures by 
an aerodynamic model that captures the salient features 
of the flow, and is valid over a large Mach number range. 
To be useful, the model must be simple enough to be 
inverted in real-time so that the airdata parameters can be 
extracted. To solve the problem of describing a complex 
flow scenario with a simple model, the FADS 
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aerodynamic model was derived as a splice of the closed- 
form potential flow solution for a blunt body, 8 applicable 
at low subsonic speeds; and the modified Newtonian 
flow model, 9 applicable at hypersonic speeds. Both 
potential flow and modified Newtonian flow describe the 
measured pressure coefficient in terms of the local 
surface incidence angle. To blend the two solutions over 
a large range of Mach numbers, a calibration parameter 
(e) was allowed for. This parameter must be empirically 
calibrated to allow for the effects of flow compression, 
body shape, and other systematic effects such as shock 
wave compression or Prandtl-Meyer expansion on the 
forebody. The resulting model is 6 

Pj= < 7 f [cos 2 (e.) + esin 2 (0 ( )] + P x (1) 

In equation 1, 0 ( is the flow incidence angle between 
the surface normal at the fth port and the velocity vector. 
The incidence angle is related to the local (or effective) 
angle of attack, ( a e ) and angle of sideslip, ( ) by 6 

cos(9 i )= cos(a f )cos(P f )cos(^) 

+ sin(P e ) sin (<)),) sin (\) (2) 

+ sin(a f )cos(P tf )cos(<(>.)sin(X, J .) 

In equation 2, the cone angle (X) is the total angle the 
normal to the surface makes with respect to the 
longitudinal axis of the nosecap. The clock angle (<j>) is 
the clockwise angle looking aft around the axis of 
symmetry starting at the bottom of the fuselage. These 
coordinate angles are depicted in figure 2. The remaining 
parameters in equation 1 are impact pressure ( q c ) and the 
free stream static pressure (P w ). Using these four basic 
parameters (a f , , q a and P^) most other airdata 
quantities of interest may be directly calculated. 

In addition to the calibration for € , the local flow 
incidence angles (a e and must be related by 
calibrations to the true free stream flow angles. These 
flow-angle calibrations account for such additional 
systematic effects as bow shock flow deflection and 
body-induced upwash and sidewash. For the X-33 a 
preliminary set of calibrations is derived using wind 
tunnel data. The wind tunnel tests and the results of the 
calibrations are discussed in detail in the Calibration of 
the Aerodynamic Model section. 

The X-33 Flush Airdata Sensing Pressure Matrix 

Since the simple model of equation 1 is derived from 
potential and Newtonian flow around a blunt body, it is 
most valid near the vehicle stagnation point. Thus, the 
most desirable location for the FADS pressure matrix is 
the reinforced carbon-carbon (RCC) nosecap which is 


used for the vehicle thermal protection system (TPS) and 
is in the stagnation region. The nosecap is hemispherical 
with a radius of 47.58 in. and extends longitudinally aft 
from the nose tip 20.8 in. Aft of the RCC nosecap an 
RCC skirt is used to protect the windward side of the 
vehicle, and a honeycomb metallic TPS is used for the 
leeward sides. 

Measurement Locations 

The number of measurements in the X-33 pressure 
matrix was selected as a compromise between the need 
to accurately measure the flow conditions at the nose and 
the cost of locating ports on the vehicle. Since there are 
four airdata states and a calibration parameter to be 
estimated, at least five independent pressure 
measurements must be available to derive the entire 
airdata state. Using five sensors to estimate the airdata is 
equivalent to a higher order spline fit and results in an 
estimating algorithm which is sensitive to noise in the 
measured pressures. Providing an additional sixth 
sensing location mitigates the noise sensitivity, increases 
redundancy options, and results in a system which gives 
overall superior performance. 

Figure 3 shows the FADS measurement locations. Five 
measurement locations are on the nosecap, and a sixth 
measurement location is on the carbon-carbon skirt. The 
layout along meridian lines allows the calculations for 
angle of attack to be decoupled from the calculations for 
angle of sideslip using the triples algorithm, the real-time 
airdata algorithm for the X-33. The mathematical 
properties of the triples algorithm will be discussed in 
detail in the Flush Airdata Sensing Estimation Algorithm 
section, and the algorithm equations are derived in the 
appendix. The system layout is designed to give good 
sensitivity for local angles of attack varying from -20° to 
45°, and angles of sideslip of up to ±20°. The nominal 
clock and cone angles of the X-33 FADS ports are 
tabulated in Table 1 . 


Table 1. X-33 FADS ports, clock and 
cone angles. 


Port no. 

<t>i , deg 

\ . deg 

1 

180 

20 

2 

270 

20 

3 

0 

0 

4 

90 

20 

5 

0 

20 

6 

0 

45 
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Pneumatic Layout 

Each X-33 flight-critical measurement subsystem 
must have a fail-operational (FailOp) capability. That is, 
the subsystem can tolerate one failure anywhere in the 
system software or hardware and still produce a usable 
result. The FADS design exploits the built-in redundancy 
of the pressure port matrix to achieve FailOp capability 
with dual redundant system hardware. This dual 
redundancy is achieved at each surface measurement 
location by installing a plug with two surface ports. 

Figure 4 shows an exploded view of the plug design. 
To survive the peak stagnation temperature of 
approximately 2000° F, the FADS plugs are fabricated 
from C-103 niobium alloy. 10 The C-103 alloy has a 
melting point in excess of 4000° F and a maximum 
working temperature of approximately 2500° F. An 
additional advantage of using niobium alloy for the 
surface plug is that the alloy has a coefficient of thermal 
expansion (CTE) which is very close to the CTE of 
carbon-carbon. 11 This match of the thermal expansion 
properties ensures that hot gas leakage around the plug is 
unlikely to occur. The relatively close match of the two 
materials expansion coefficients also ensures that undue 
thermal stresses are not imparted to the RCC nosecap 
when high heating loads are applied. The plug design 
was derived from an adaptation of the Space Shuttle 
Entry Airdata System 3 which used niobium inserts with 
excellent success. 

On each plug the two surface ports are connected to 
the pneumatic tubing using a high-temperature titanium 
alloy braze. The individual surface ports are plumbed to 
the measurement sensors using approximately 8-12 ft of 
pneumatic tubing with an approximate inside diameter 
of 0.21 in. The pneumatic response of this arrangement 
was analyzed extensively along the X-33 flight 
trajectory. Based on these analyses, it was determined 
that the pneumatic latencies in the system would not 
introduce significant errors into the airdata estimates. 

Computational Architecture and Redundancy 
Management 

A schematic of the computational architecture and 
redundancy management scheme for the X-33 FADS is 
presented in figure 5. The design of the X-33 FADS 
system provides a total of twelve surface pressure 
measurements; however, the pressures from the dual 
redundant pressure ports are always analyzed 
independently. Defining data flow path I as the set of 
computations which use the grouping of the six upper 
and outboard pressure measurements on each plug, and 


dataflow path II as the set of computations which uses 
the six lower and inboard pressures, then the FailOp 
capability of the system will always be ensured by 
selecting the computational path with the minimum 
mean square fit error. The mean square fit error is 
computed as 

x p.) 

i = \ 

^2 + esin 2 0.(a e , $ e )] + A*,] 2 

X* 

i- I 

where, the A designates a quantity estimated by the 
FADS algorithm. This redundancy management scheme 
selects the system with the best overall fit consistency, 
and allows for a soft sensor failure (one which is not 
detected by the hardware diagnostics) to occur without 
degrading the performance of the system. 

If the mean square error (MSE) of the output flow path 

is normalized by an expected population variance (that 

is, by the expected range of fit error that is allowable for 

a system with no failures) then the MSE becomes 

distributed as y} and is a good indicator of the absolute 

system health. This chi-squared approach was developed 

J p n 

and successfully employed by Whitmore, et al * on 

previous flight tests. Expected values for the fit error 

ranges are still to be determined from additional wind 

tunnel tests. 

The Flush Airdata Sensing Estimating Algorithm 

The aerodynamic model of equations 1 and 2 is 
inherently nonlinear in terms of the state parameters and 
can not be directly inverted. A nonlinear regression 
method in which the equations are recursively linearized 
and inverted using iterative least squares has been 
developed and successfully flight tested. 14 Problems 
with the stability of the algorithm were encountered, 
however, and the algorithm required special software 
patches to maintain stability. Because the FADS is to be 
used for closed-loop flight control on the X-33, the 
nonlinear regression algorithm w ; as determined to be too 
risky and was abandoned. 

A new solution algorithm was developed for the X-33 
FADS. A better solution algorithm is offered by taking 
strategic combinations of three sensor readings that will 
decouple the flow angularity states from the static and 
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impact pressure states and the aerodynamic calibration. 
Detailed derivations of the FADS estimator equations are 
presented in the appendix. For completeness, portions 
necessary to explain the flow of the algorithm are 
restated here. 

Angle of Attack Triples Estimator 

By strategically taking the differences of three surface 
pressures (referred to in this paper as “triples”) the 
pressure related states ( q c , and e) are eliminated 

from equation 1 . The resulting pressure equation is 

r, * cos 2 Qj + Tji cos \ + T k j cos 2 0 ( = 0 (4a) 

where 

r ,A = (Pi ~ Pk )' r ji = ( Pj ~ Pi >- r *, = - P? (4b) 

and pi , Pj, and p k are the pressures used in the data triple. 
In equation 4 the local angle of attack can be decoupled 
from the local angle of sideslip using only pressures 
aligned along a vertical meridian where 0 = 0° or 180°. 
In this geometry arrangement, terms related to angle of 
sideslip fall out of equation 4. For la I < 45° the solution 
for a is 


a = 




(5a) 


where 

/l = r, * sin + r^sin \. + sin \ 

B = r ik costyjSinXjCoshj 
+ T^cos^sin Xf.cosXj, 

+ r^cosfysin/^cos^j 


(5b) 


When la! > 45° the correct solution is given by the 
complement of equation 5(a) 


a = 




(5c) 


clock and cone angles corresponding to these triples 
combinations are presented in Table 2. 


Table 2. Angle of attack triples, clock and cone angles. 


Triple no. 

4> i- 
deg 

deg 

deg 

b 

deg 

<t >*’ 

deg 

k*' 

deg 

1 

0 

0 

0 

20 

180 

20 

2 

0 

0 

0 

20 

0 

45 

3 

0 

0 

180 

20 

0 

45 

4 

0 

20 

180 

20 

0 

45 


The output angle of attack estimate is determined as 
the mean of the values computed using the four 
individual triples. This averaging procedure provides a 
measure of noise rejection for the estimator. Clearly, if 
one of the ports along the vertical meridian is deemed 
unusable and is weighted out of the algorithm, then only 
one valid triple remains for computing the angle of 
attack. 

Angle of Sideslip Estimator 

Once the local angle of attack has been estimated, then 
the angle of sideslip may be evaluated using any 
combination of the available ports, other than the 
obvious set in which all three ports lie on the vertical 
meridian. The result, derived in the appendix, is a 
quadratic equation in tan P . 


Atan 2 p + 2Z?'tan P + C = 0 


(6a) 


where 


/t’-iratjUr^+r^V 


B' = \ r ik ajbj + T }i a k b k + T kj a i b l 


(6b) 


c = \ r ik a j 2+r ji a k 2+r kj a , 2 


For the X-33 there are four possible independent 
combinations of a-triples on the vertical meridian. The 


and 
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a {ijk} = cosacosA {jW 

+ sinasin^ /;jt jCos<t) { ^} (6c) 

b {ijk] = smX {ijk} sin ${ijk} 

Equation 6 has two solutions. However, unlike the 
solution for angle of attack these solutions cannot be 
reduced to a single obvious choice for the angle of 
sideslip roots. Determining which root is correct 
depends on the port arrangement used to determine the 
angle of sideslip. Since only four angle-of-sideslip ports 
are required for a redundant measurement system, the 
number of possible angle-of-sideslip triples to be 
considered is reduced to a more manageable number by 
removing the upper 20° port (port no. 1 ) on the vertical 
meridian. When this port is eliminated the number of 
available triples reduces to seven. The list of triples to be 
considered is tabulated in Table 3. 


Table 3. Angle of sideslip triples, clock and cone 
angles. 


Triple 

no. 

4>r 

deg 

deg 

*>• 

deg 

deg 

deg 

K 

deg 

1 

0 

0 

90 

20 

270 

20 

2 

0 

20 

90 

20 

270 

20 

3 

0 

0 

90 

20 

0 

20 

4 

0 

0 

270 

20 

0 

20 

5 

0 

45 

90 

20 

270 

20 

6 

0 

0 

90 

20 

0 

45 

7 

0 

0 

270 

20 

0 

45 


{ A., <j) } = {0°, 45°} is used only when a e is outside of 
the range from 17° to 20°. Similarly, only when a e is 
between the range from 17° to 20°, then is port number 5 
{>.,<(>} = {0°, 20°} used as the fourth angle-of-sideslip 
port. 

As with the angle-of-attack algorithm, the output 
angle of sideslip estimate is determined as the mean of 
the values computed using the four individual triples. 
Clearly, if one of the ports along the vertical meridian is 
deemed unusable and is weighted out of the algorithm, 
then only one valid triple remains for computing the 
angle of sideslip. 

Mach Number. Static Pressure, and Impact Pressure.. 

Es timato r 

Once the values of a e and (J f have been determined, 
then the incidence angles at all of the ports can be 
evaluated, and only e , , and q c remain as unknowns 

in the pressure equations. As derived in the appendix, e 
is implicitly a function of P ^ and q c and the resulting 
system of equations is nonlinear. The solutions for 
and q c must be extracted iteratively. Defining the 
matrices 


f cos 2 (0j) > 
e <J ' ) sin 2 (0 1 ) y 


<7i 


. . . 0 


M 


U) 


, Q = 


(7) 


( cos 2 (0 6 ) N 

V+ e u, sin“(0 6 ) J 


0 . . . q 6 


As developed in the appendix, port arrangements 2 
and 5 have mathematical singularities at 0° for certain 
values of a e . Near these singularities all three pressures 
in the triple have nearly the same value, and the triples 
equations become ill-conditioned. Away from these 
singularities, the solution for is given by the root 
whose absolute value is closest to zero. However, as 
derived in the appendix, near the singularity the solution 
is indeterminate, and the proper solution cannot be 
distinguished. 


the original flow model (eq. 1) can be used to develop an 
iterative estimator of the form 


<*c 

i A 


°jo + /> 




-l 




P i 


P 6 


(8) 


For port arrangement 2 the singularity occurs at 
0° a e . For port arrangement 5 the singularity is at 
18.2° a e . For the X-33 FADS design, which requires 
four ports for the estimator, the singularity problem is 
managed by a simple algorithm in which port number 6 


The subscript (j) refers the result of the / th iteration. 
The notation M {J) refers to the matrix of equation 8, with 
e being evaluated, using the values for Mach number 
resulting from the previous iteration. The q x terms are 
weights which have a nominal value of 1.0. Setting the 
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value of q, to zero, weights the i'th pressure reading out 
of the algorithm. Equation 8 is derived in detail in the 
appendix. 

Given q c and Mach number can be computed 
using normal one-dimensional fluid mechanics 
relationships. Subsonically, Mach number can be 
calculated directly using isentropic flow laws, where for 
y = 1.4 


£ = [1 + -2mL] 3 ' 5 - 1 (9a) 


Supersonically, the solution is computed using the 
Rayleigh pitot equation which is derived from adiabatic 
normal shock wave relationships. 15 For y = 1.4 


Qc 


P 


oo 


166.92mL 

2 ^5 

[7Mi-l] 


(9b) 


Equation 9(b) is solved using a Taylor series 
expansion and a reversion of series to solve for Mach 
number. 16 


For high Mach numbers, a very significant 
temperature rise occurs across the bow shock wave and 
the fluid no longer behaves as a calorically perfect gas. 
The large temperature rise across the shock wave causes 
vibrational modes of the air molecules to excite, which 
draws energy away from the flow. Thus the stagnation 
temperatures actually encountered are considerably 
lower than those which would be computed if perfect gas 
calculations with y — 1 .4 were used. 17 


As a result of this nonadiabatic flow behavior, some 
error is expected to occur in the Mach number 
calculation using equation 9(b). Fortunately, because 
pressure is a mechanical quantity and depends primarily 
upon the mechanical aspects of the flow, influences 
caused by high-temperature gas properties are 
secondary. Numerical analyses of equation 9(b) have 
shown that the error introduced is less than 0.2 percent 
for Mach numbers below 4.O. 15 ’ 18, 19 This error is 
acceptable for the X-33 design. 


analyze the stability of nonlinear equations with two 
unknowns; however, a linearized stability analysis 
determines the behavior of the system with respect to 
small disturbances. For a discrete iteration, eigenvalues 
with magnitudes less than 1.0 indicate stability 20 for the 
corresponding eigenmode. As the magnitudes of the 
eigenvalues approach unity, the system is relatively less 
stable, and more iterations are required for convergence. 
Eigenvalues with magnitudes greater than unity indicate 
that the iteration will diverge. 

For the FADS algorithm there are two eigenvalues. As 
derived in the appendix, these eigenvalues are primarily 
a function of Mach number and, to a much lesser extent, 
angle of attack. The individual values of q c and did 
not affect the stability of the iteration, only their ratio is 
significant. This conclusion is very important, as it 
allows the stability characteristics of the system iteration 
to be analyzed independently of specific trajectory. 
Figure 6 shows the eigenvalue magnitudes as a function 
of Mach number and angle of attack. One eigenvalue is 
always zero and has no effect on the stability of the 
system. The other eigenvalue has a nonzero magnitude 
and determines the convergence properties of the 
algorithm. 

Subsonically, the nonzero eigenvalue is stable, but 
approaches neutral stability at Mach 1.0. This nearly 
neutral stability causes the algorithm to take a larger 
number of cycles to converge transonically. As presented 
in the Evaluation of the System Accuracy section, this 
slow convergence also causes the overall system 
accuracy to be somewhat degraded near Mach 1.0. For 
moderate supersonic Mach numbers, the eigenvalue 
magnitude drops rapidly and the algorithm again 
becomes stable. As the Mach number increases beyond 
Mach 4, the stability of the algorithm is reduced until 
eventually the iteration becomes unstable beyond 
Mach 8. The likely cause of this reduced stability is that 
equation 9(b) has become so ill-conditioned that small 
changes in the ratio of q c to P ^ produce large changes in 
the free stream Mach number. Fortunately, below 
Mach 4, where airdata are required by the vehicle flight 
control system, iteration stability does not appear to be a 
problem. 


Numerical Stability of the Flush Airdata Sensing s 
Algorithm Iteration 

Since equations 8 through 9 are to be implemented as 
part of a real-time airdata estimation algorithm, it is 
essential that potential instabilities in the iteration 
scheme be identified. In general it is not possible to 


Al gorithm Startup 

At very high Mach numbers, even though the and 
solutions will still be valid, the stability problems 
presented in the last section require the iterative 
algorithm to be shutoff to avoid numerical divergence. 
When this algorithm shutoff is applied, a reliable method 
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must be available for bringing the algorithm back on-line 
when the Mach number decreases back to safe limits. 
One approach relies on the use of inertial Mach number 
to initialize the algorithm. Another simple, but effective 
approach requires no outside information for 
initialization and is presented here. The method 
presented here was developed and flight tested by 
Whitmore, et al 6 (1995). 

A block diagram of the startup algorithm is depicted in 
figure 7. A prescribed startup Mach number is 
preselected and hard-coded. For the X-33 trajectory the 
default startup Mach number is 4.0. When the startup 
algorithm is initiated, the transducers are polled and an 
initial pressure data set (the target) is obtained. Using the 
target pressure set, a e and are evaluated 

independently of Mach number using the triples 
algorithm. With the resulting values for a e and and 
the target pressure data, the FADS algorithm iteration is 
used to compute a new estimate for static pressure, Mach 
number, impact pressure, and e. 

The aerodynamic model is then used to generate a 
predicted pressure distribution corresponding to the 
computed airdata set. The resulting biased pressure 
distribution is linearly perturbed toward the target set by 
a small increment (currently 1/100 of the distance 
between the target set and the predicted pressure set). 
Using the perturbed value, a new estimate of the Mach 
number is evaluated using the FADS iterative algorithm. 
A new pressure prediction corresponding to the 
perturbed airdata state is evaluated and the progression 
along the line from the predicted pressures to the target 
pressures is repeated again until the target set is reached 
or the maximum number of increments have been 
performed. 

Calibration of the Aerodynamic Model 

There are three calibration parameters which must be 

evaluated for the X-33 FADS system: the position error 

(e), the angle of attack flow correction angle (5a), and 

the angle of sideslip flow correction angle (5(3). For the 

preliminary design these calibration data were obtained 

from wind tunnel data. A two-percent model of the full- 

scale X-33 was instrumented with pressure taps 

corresponding to the desired X-33 locations and then 

tested over a wide Mach number, angle of attack, and 

angle of sideslip range. To obtain the preliminary 

calibrations, the model was tested in two NASA Langley 

1 

wind tunnels (1) the 16-Foot wind tunnel and (2) the 
Unitary Plan wind tunnel (UPWT). 22 A description of 
the test model, the test instrumentation, and the test 
conditions follows. 


X-33 Wind Tunnel Model and Test Instrumentation 

The two-percent X-33 model was manufactured for 
measuring total vehicle forces and moments and 
measuring nosecap surface pressures. The 21 nosecap 
surface pressures included the locations of the X-33 
FADS ports and were used to determine the FADS 
aerodynamic calibrations. The wind tunnel model was 
machined from a solid piece of aluminum and the nose 
cone was made of stainless steel to minimize damage 
from tunnel contaminants. All FADS pressure ports were 
drilled normal to the surface with a diameter of 0.04 in. 

The surface pressures were sensed through an 
electronically scanned pressure (ESP) module which 
produces a time-multiplexed analog output with up to 
32 channels available for each module. The ESP was 
located approximately 5 in. aft of the nosecap and was 
plumbed to the various surface ports using flexible 
pneumatic tubing. Because of the short line lengths 
involved, pneumatic lag in the pressure tubing was 
considered negligible. The ESP sensor was a ±10 lb/in“ 
differential module with a manufacturers accuracy of 
better than ±0.1 percent of the full-scale reading. The 
time-multiplexed analog outputs from the ESP were 
tagged and sampled using a 16-bit analog-to-digital 
(A/D) conversion system. 

During the wind tunnel tests, zero-shifts in the ESP 
module were frequently monitored and later applied to 
the data. Additionally, spare ESP ports were plumbed 
together to monitor the drift in the individual sensors in 
the ESP unit during the tests. If, for a particular test 
point, the drift rates were observed to be excessive, then 
the test point was repeated. Test points in which any ESP 
drift exceeded 2.0 lb/fr were not used. The reference 
pressure source for ESP was located outside of the wind 
tunnel and was sensed using a highly accurate absolute 
pressure transducer. The reference pressure transducer 
featured a capacitive transduction technology, with an 
accuracy of ±0.25 lb/ft 2 . Outputs from the sensor were 
sampled using the 16-bit A/D system. 

Wind Tunnel Test Facilities and Test Conditions 

The Langley 16-foot wind tunnel is a closed-circuit, 
single-return, continuous-flow', atmospheric tunnel with 
a Mach number range from 0.2 to 1.3. The tunnel has an 
octagonal, 16 ft slotted test section. Run limitations are 
often imposed near transonic Mach numbers because of 
the lack of porous walls that absorb reflected shock 
waves. The small size of the X-33 model, however, 
allowed Mach numbers of 0.95 and 1 .05 without 
reflected shock wave interference. Angles of attack 
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( 10 ) 


ranged from -10° to 24° and angles of sideslip from -10° 
to 16°. Angle of sideslip sweeps were done at 0° and 8° 
angle of attack, and angle of attack sweeps were done at 
-8°, -5°, 0°, 5°, and 8° angles of sideslip. 

The UPWT is a closed-circuit, continuous-flow, 
variable density supersonic wind tunnel with two 4 x 4 ft 
test sections. One test section has a Mach number range 
of 1.5 to 2.9 and the other has a range of 2.3 
to approximately 4.5. The angle of attack ranged from 
-10° to 26° and angle of sideslip from ±10°. Angle-of- 
sideslip sweeps were done at 0° and 8° angle of attack. 
Angle of attack sweeps were done at -8°, -5°, 0°, 5°, 
and 8° angles of sideslip. The Reynolds number was 
nearly constant at 2 x 10 6 perft. As with the previous 
16-ft tunnel tests, runs were conducted with and without 
a grit ring installed aft of the FADS ports, with little 
variation in the pressure distribution observed. 

Table 4 summarizes the ranges of conditions and 
configurations which were tested. 


5a = A Q [MJ+A { [MJa e 

+ A 2 [MJaJ + A 3 [MJaJ 

The true angle of attack is evaluated by subtracting the 
correction factor (5a) from the indicated value. The 
calibration data and the corresponding curve fits for 5a 
are presented as a function of Mach number in figure 8. 
The polynomial fits to the wind tunnel data are excellent. 
Also notice that the curves flatten out with increasing 
Mach number, a result of diminishing flow deflections at 
increasing Mach numbers. 

Angle-of-Sideslip Calibration 

In a manner identical to the angle-of-attack 
calibration, the angle-of-sideslip calibration relates the 
correction to the local angle of sideslip sensed by the 
FADS to the free stream value. The wind tunnel residual 
data were fit with third order polynomials in P e with the 
coefficients scheduled as a function of Mach number. 
The resulting calibration function is 


Table 4. Range of wind tunnel test conditions. 


Tunnel 

Mach range 

a, 

deg 

p. 

deg 


0.25 to 

-10 to 

-10 to 

16-ft 

1.2 

24 

16 


1 .6 to 

-10 to 

-10 to 

UPWT 

4.50 

26 

10 


Evaluation of the FADS Calibration 

The calibration parameters were estimated from the 
wind tunnel data using measured pressures to estimate 
a e and $ e by means of the triples algorithm. Results for 
a € and $ e were used with wind tunnel reference 
conditions for M^ and q c to predict surface pressures at 
the measurement locations. The residuals between the 
predicted and measured pressures were then used to 
calculate e for each test point using linear regression. 
Finally, trends in the residuals were curve fit to produce 
the calibrations. A description of the three component 
parts of the calibration follows, in detail. 

Angle-of- Attack Calibration 

The angle-of-attack calibration relates the local angle 
of attack to the free stream value. Residuals between the 
wind tunnel reference conditions and the estimates for 
a € were curve fit with third order polynomials in a e to 
give the correction factor (5a). The coefficients of the fit 
were scheduled as a function of M^. The resulting 
calibration function is 


5P = + 

+ b 2 [mjp* + b 3 [ m jpJ 


( 11 ) 


The calibration data and the corresponding curve fits 
for 6(3 are presented in figure 9. 

Position Error (e) Calibration 

The position error calibration parameter adjusts the 
FADS pressure model for changes in the incidence angle 
and Mach number. The wind tunnel data were curve fit 
with second order polynomials in angle of attack and 
angle of sideslip, and the coefficients scheduled as a 
function of Mach number. The resulting calibration 
model is 


E = e M [HJ + £ a JMJa e + e a JMJaJ 
+ ep l [MjP e + e Pi [M M ]pJ 


( 12 ) 


The calibration curves for £ at zero p e are presented 
in figure 10(a). The variations of e with $ e at zero a e 
are shown in figure 10(b). The parameter e M (at zero a e 
and P ) is plotted as a function of the free stream Mach 
number in figure 10(c). 

It is interesting to note that the trends shown in 
figures 10(a) and 10(b) are remarkably similar. In both 
figures, 10(a) and 10(b), the data form a family of 
parabolic curves. The peak of each parabola (e M ) 
increases in magnitude and shifts to increased angle of 
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attack as Mach number increases. The curves in figures 
10(a) and 10(b) have negative concavity subsonically, 
almost no concavity transonically, and positive 
concavity supersonically. Also, at close to supersonic 
Mach numbers, e M approaches zero, an observation 
which confirms that the FADS flow model is consistent 
with modified Newtonian flow theory. 

Evaluation of the System Accuracy 

The accuracy of the calibrated FADS system was 
analyzed using an analysis-of-variance 23 method. In this 
method, the calibrations for e, 5a, and 8(3 and also the 
measured pressure distributions were used to generate 
airdata estimates with the FADS algorithm for each wind 
tunnel test point. Assuming that the measured wind 
tunnel reference conditions represent a truth set, 
quantitative measures of the overall system accuracy are 
a product of the residuals between FADS algorithm 
estimates and the wind tunnel reference conditions. The 
residuals include errors in the wind tunnel reference 
conditions, and are considered to be representative 
measures of the absolute system accuracy. 

Note that the calibration curves were derived using 
only data from angle of attack sweeps at zero angle of 
sideslip and angle of sideslip sweeps at zero angle of 
attack. The error analysis, however, includes residual 
data that were obtained during test points which 
included both nonzero angle of attack and nonzero angle 
of sideslip. 

Figure 1 1 shows the available wind tunnel database 
for angle of attack and Mach number and typical 
Michael Army Air Field (MAAF) and Malmstrom Air 
Force Base (MAFB) flight trajectories. Clearly, many of 
the wind tunnel points are far outside of the actual flight 
conditions to be encountered during the X-33 flight tests. 
Introducing these unrealistic data into the residual 
analysis would likely give error estimates which are too 
large. Thus, the following procedure was used to reduce 
the wind tunnel data sets to a more realistic and 
manageable size. 

1. The outer boundaries of the minimum and 
maximum angles of attack and angles of sideslip to 
be encountered as a function of Mach number were 
determined for typical MAAF/MAFB launch and 
TAEM flight trajectories. 

2. Next, the union of the outer boundaries was 
determined and ±5° angle of attack and ±3° angle 
of sideslip were added to give the domain space 


over which the residuals would be analyzed. This 
domain space allows for considerable variation 
away from the nominal trajectory. 

y 

3. All test points which exhibited more that 2.0 lb/ft 
of drift in the ESP readings were eliminated from 
the domain space. 

4. Within this domain space, root-mean-square 
(RMS) averages of the residuals for the angle-of- 
attack and angle-of-sideslip test points were 
evaluated as a function of Mach number. 

5. Finally, the resulting RMS averages were fared and 
tabulated as a function of Mach number to produce 
a set of three simple 1 -G error models with Mach 
number as the independent variable, and SM, 8a, 
and 8(3 outputs. Figure 12 displays these RMS 
error estimates as a function of Mach number. 

Figure 12(a) shows the percentage of Mach number 
error and the accuracy requirement for the X-33 FADS. 
The Mach error shows a sizable error jump at transonic 
speeds, but does not exceed the error limit. The angle-of- 
attack and angle-of-sideslip error plots are displayed 
with the accuracy requirements in figures 12(b) and 
12(c). Below Mach 4.0 the errors for both parameters are 
within the accuracy requirements of the X-33 airdata 
system. 

Note that in examining figure 12(a) in more detail, 
growth in the Mach number error percentage (both 
transonically and above Mach 2.5) corresponds closely 
to the growth in the nonzero eigenvalue (fig. 6). This 
result is not entirely surprising. The eigenvalues 
represent a measure of the numerical conditioning of the 
system. Transonically this numerical ill-conditioning is 
caused by the inflection point at which the flow sw itches 
from subsonic to supersonic. This inflection point 
corresponds to a point of maximum slope in the e M 
Mach number curve (fig. 10(c)). At supersonic speeds 
this ill-conditioning is caused by the rapid growth in the 
strength of the bow shock as Mach number increases. 
The result is that large changes in the free stream Mach 
number produce only small changes in the local pressure 
distribution. Conversely, small changes in the measured 
pressure distribution make for large changes in the 
estimated free stream Mach number. This ill- 
conditioning with respect to Mach number is an 
unavoidable consequence of the flow physics and is 
common to all airdata measurement systems at transonic 
and supersonic Mach numbers. 
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Summary and Concluding Remarks 

The X-33 FADS design is presented in this paper. The 
FADS concept, in which airdata are inferred from 
nonintrusive surface pressure measurements, does not 
require stagnation of the surface flow. This innovation 
extends the useful range of the airdata measurement 
system to the hypersonic flow regime. The FADS system 
uses a matrix of pressure orifices on the vehicle nose to 
estimate airdata parameters. 

The FADS system design exploits the built-in 
redundancy in the pressure port matrix to achieve FailOp 
capability with only dual redundancy in the system 
hardware. The design produces two independent 
measurement paths for the airdata and the redundancy 
management scheme selects the measurement path with 
the best overall fit consistency. This scheme allows for 
failure, either detected or undetected, to occur anywhere 
in the system without degrading the overall performance. 

Presented here is the aerodynamic model that relates 
the airdata parameters to the measured pressures. Also 
presented is the derivation of the estimating algorithm. 
The mathematical properties of the algorithm were 
analyzed in detail. At near hypersonic Mach numbers, 
the stability of the algorithm is reduced because the bow 
shock wave has become so strong that small disturbances 
are highly amplified. Eventually the iteration becomes 
unstable, at close to Mach 8 and beyond. This ill- 
conditioning with respect to Mach number is an 
unavoidable consequence of the flow physics. 
Fortunately, below- Mach 4, where airdata are required 
for the vehicle flight control system, potential 
mathematical singularities and iterative instabilities in 
the nonlinear estimating algorithm can be prevented by 
using simple logic. 

Because of the problems with stability at high 
supersonic Mach numbers, the iterative algorithm must 
be shutoff in order to avoid numerical divergence 
problems. A reliable method for bringing the algorithm 
back on-line was presented. This simple creeping 
solution method insures that the algorithm does not 
exceed the small perturbation limits of the stability 
analysis. 

There are three calibration parameters which must be 
evaluated for the X-33 FADS system: the position error 
(£), the angle of attack flow correction (5a), and the 
angle of sideslip flow correction (5P). For the 
preliminary X-33 design these data were evaluated from 
wind tunnel calibrations. A two-percent full-scale X-33 


model was instrumented and tested over a variety of 
angles of attack and Mach number ranges in two separate 
wind tunnels. Data were obtained for Mach numbers 
ranging from 0.25 to 4.5, angles of attack from -10° to 
25°, and angles of sideslip from ±10° to 16°. 

The accuracy of the system was evaluated by applying 
the calibrations and the FADS algorithm to generate 
airdata estimates for each wind tunnel test point using 
the measured pressure distribution. Assuming that the 
wind tunnel reference conditions represented a truth set, 
the statistical accuracy was evaluated by computing 
residuals between the reference set and FADS estimates 
and then compiling the results into a set of 1-0 error 
models. The error analysis concluded that the FADS 
meets the accuracy requirements for the X-33 vehicle. 

Appendix: Mathematical Analysis of the 
Flush Airdata Sensing Algorithm 

As discussed in the main text, the airdata states are 
related to the surface pressures by an aerodynamic 
model that captures salient features of the flow, and is 
valid over a very large Mach number range. The 
resulting model takes the form 

Pi = < 7 c .[cos 2 ( 6 ( ) + esin 2 (0 ( )] + P m (A-l) 

where 

cos(0 ( ) = cos(a f ,)cos(P f ,)cos(^ / ) 

+ sin(P f )sin(<j> ( )sin(>. ( ) (A-2) 

+ sin(a <1 )cos(P ( ,)cos(())-)sm(^ 1 ) 

Defining the terms 

a = cosa cosk, + sina sin^costi 
' ' (A-3) 

b t - sin^sinfy 

equation A-2 reduces to 

cos(G-) = a-cos(P^) + bjSinCP^) (A-4) 


Using equation A-l to take strategic combinations of 
three pressures (triples) the impact pressure, static 
pressure, and calibration parameter are decoupled from 
the local angle of attack and angle of sideslip. 
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Pi-Pj 


Pj-Pk 


J _ 


qA cos 2 0 ( + Esin‘0 | |- + P o , 


2 2 

q< cos 0 ; + esin 0,1- + P„ 


<7 C <! cos 2 0y + esin 2 0 y + P a 


2 2 

< 7 C <| cos 0 A . + esin 0 k \ + P a 


(A-5) 


j(l -e)cos 2 0 ( + e 


(1 - e)cos 2 0 ; + e| 

|(1 -e)cos 2 0 ; + ej 

H 

|( 1 - z)cos 2 & k + e j 


2 2 

COS 0 . - cos 0 ^ 

“ 2 2 

COS - COS 0 £ 


Expanding equation A-5 and defining 

T ik s Pi-p* r ji-Pj-Pi’ r kj s Pk~Pj < A -o 


2 

r^[cosacos^ + sinasin^cost^] 

2 

-f iy cosacos^ + sinasinX A ,cos<j>j.] (A-9) 

2 

+ ry cosacosA.; + sinasin^cos<t> ( r = 0 
Factoring cos a (a*± 90°) out of equation A-9, 
r ik [cosXj+ tan a sin ^ cos (J)^] 2 

2 

+ rycosA.** tanasin A^.cos<j)J (A- 10) 

2 

+ ry cosA* f + tanasin A^cos<}) ( ] = 0 

Expanding the squares, the result is a quadratic 
expression in tana, 

2 2 2 2 

[T^sin Xj + iy sin A^ + r^sin X y ]tan a 

+ 2 [ cos Xj sin X- cos + r^cosA^sinA^cos^ 

+ r^cosA^sinA^cos^ ]tana ^ ^ 

+ [r jk cos 2 Xj + T ji cos 2 X k + r k jCOS 2 X i ] = 0 

Equation A- 11 can be reduced further by noting that 

2 2 
cos X = 1 - sin X and 


equation A-5 reduces to 

r,*cos 2 e ; + r y/ cos\ + ryos 2 e. = o (A-7) 

Substituting in from equation A-4, dividing by 
2 

cos P (P *±90°) and collecting terms, equation A-7 
becomes 

r *[a + ft tan P] 2 + Y [a k + ft^tanP] 2 

J 1 1 (A-8) 

+ r kj [a l + b l tanP] 2 = 0 

In equation A-8, the local angle of attack and angle of 
sideslip have been written as a function of the measured 
pressure only, the impact pressure, £, and static pressure 
have been decoupled. 

Angle of Attack Solution 

Angle of attack is decoupled from angle of sideslip 
using a meridian solution, in which, along the vertical 
meridian <|) = 0, ±180°, and equation A-8 reduces to 


+ r ji COs2l k + V kj cos \ 

= r ik (]-sin 2 \ j ) + r ji (\-sm\) 

+ r,/i-sin%) = r lk + r jl + r kj 

- (T lk sin 2 Xj + r ; ,sin\ + T kj sm\) < A ' 12) 

= (Pj-P k ) + (Pj-Pi) + (P k -Pj) 

- (Y jk sin 2 X j + r ji sm\ + Y kj s\n\) 

= -( r ,i.sin 2 ^ + r j( sin 2 ^, + r^sin%) 

Substituting equation A- 1 2 into equation A- 1 1 , 

[r iJt sin\ y + r^sin 2 ^ + r^sin 2 X ( ][tan‘a - 1] 

+ 2[r ( . t cos^sinX ; .cos()) 7 (A ., 3) 

+ TycosA^sinA^cos^ 

+ r^cos^sinA-cos^Jtana = 0 
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Defining 


Defining 


2 2 2 
A = (r.^sin Xj + T ; ..sin X^ + r^sin X ) 

B - (T ik cos$jSinXjCosXj 

+ r^cos^.sinX^cosX^ + r^cos()).sinX.cosX.) 

Equation A- 14 is written in the simple form 


(A- 14) 


A[tan 2 a- 1] + 2fltana = 0 


(A- 15) 


The correct solution is extracted from equation A- 15 
by noting that 


. 2 2 

2 , sin a - cos a 

tan a - 1 = 

cos a 


1 - 2 cos a 
2 

cos a 


cos 2a 
cos a 


(A- 16) 


Substituting equation A-^6 into equation A- 15, and 
multiplying through by cos “a , 


A’ = (T ik bj 2 +r ji b k 2 + r kj b i 2 ) 

B' = {T ik a ] b ] + T Jt a k b k + T^afo) (A-21) 
C = (T ik a 2 +V jpf+T^a 2 ) 

Equation A-20 reduces to a quadratic equation in 
tan [5 

A'tan 2 ^ + 2B' tan(J + C = 0 (A-22) 

The two solutions to equation A-22 are 


tan (ip tanp 2 



(A-23) 


Root Selection Criteria for Angle of Sideslip 

The relationship of tan [ P j ] to tan[P 2 ] is developed 
by taking the product of the roots 


-/I cos 2a + 2£cosasina 
= -4cos2a + 5sin2a = 0 

and the solution for a is 


a 



~A~ 

_B_ 


(A- 17) 


(A- 18) 


For la I < 45°, equation A-18 picks the correct root; 
when la I > 45° then the correct root is given by the 
complement of equation A- 1 8 


a = 


1 

2 


n - tan 



(A- 19) 


An gle of Sideslip Solution 

Given the solution for angle of attack, a procedure 
similar to that of the previous section is followed for 
angle of sideslip. Expanding the squares in equation A-7 
and solving for like terms in tan P, 



Rearranging equation A-24 gives the relationship of 
the two roots 



For several specific port arrangements relevant to the 
X-33 design, the relationship between the two roots can 
be written in closed form. Three of these special 
arrangements will now be discussed. 

Root Selection Criteria for the Triples on Lateral 
Meridian 


\Tik b ) + r ji b l + r^J]tan 2 p 
+ 2fT >/, + 1 >*V + (A ‘ 20) 

+ r ,k a ) + r j, a l + T k/i = 0 


Along the lateral angle-of-sideslip meridian, 
the outer port clock angles are ±90°, and 
ctijf, = cosacos X ijk , b- { j k = ± sinX^, thus equation 
A-24 reduces to 
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[*]= 


C"] = T ik a j +r ji a k + r kj a i 

r;.6, 2 + r,,fc, 2 + r w 6 f 


/A"; 


jr* 


kj i 


= cos a 


£^fos ^ + r ; ., cos ^ + r^ co5 \ 
r i* sin % + r ;,sin\ + T^shi 


(A-26) 

2-v ^ 


= —cos a 


Thus, along the lateral meridian the relationship 
between the two solutions of equation A-22 is 


C 2 
= cos a 
A 


o 2' 

-T cosX + T ..(cosX, + tan a sin X.) 

J H 2 11 


2 

cos a 


2 . | C0S ^ 2 
sin X . cos a 

2 2 2 
_ cos acosX^ - cos (a-Xp 

sin X- 


-r. /S in 2 x. 


1 2 
( cos a cos X, + sin a sin X, ) 


(A- 30) 
•*> 


tan [3 j = cos a tan 




(A-27) Substituting equation A-30 into equation A-25, the 
relationship between the two roots becomes 


The solution for P 2 ls plotted as a function of P t and 
the local angle of attack in figure A- 1 . Even for angles of 
attack as high as 50° the two angles of sideslip solutions 
are nearly 70° out of phase; so for triples along the 
lateral meridian, selecting the root with Ip I closest to 
zero will always produce the correct solution for p . 


tanp, = 


cos 2 acosXy 2 - cos 2 (a - X k ) 
sin'X, 


(A-31) 


x tan 


L2 


-P 2 


Off-Meridian Triples. Symmetric Arrangement 

Assume ports / and j lie on the lateral meridian and 
port k lies on the vertical meridian, then 

<[)• = - => a- = cosacosX^ 

b { - sinX / 
tt % 

<t >j = 2 =>a j = cosacos Kj (A-28) 

bj = -sinX^ 

§ k = a k = cosacosXj. + sinasinX^ 

h= o 

C 

and for a symmetric arrangement, X- = X^, and — 
reduces to 

C 2 
— = cos a 

A (A-29) 

^ 2” 

cosX -(T ik + T k j) + r ; ,(cosX t + tanasin^.) 

sin\ j(T ik + r k j) 

Noting that r ik + T kj = - T ji =» equation A-29 can 
be reduced further 


For the X-33 design, there are two possible symmetric 
off-meridian angle of sideslip triples combinations, 


{X,, X r X k } = {20°, 20°, 45°}, and 
{X ( , X r X k ] = {20°, 20°, 20°} 


(A-32) 


For the first set of coordinates {X^X^, \ k } = 
{20°, 20°, 45°} the solution for p 2 is plotted as a 
function of Pj and the local angle of attack in 
figure A-2(a). Notice that as a approaches 1 8.207°, the 
angle of attack at which all of the incidence angles are 
equal, a singularity exists at P =0. Thus, for angles of 
attack near 18.207°, and angles of sideslip values near 
zero, the system is highly ill-conditioned. This ill- 
conditioning makes selecting the correct root 
impossible. Figure A-2(b) presents the same analysis as 
figure A-2(a), except now the second set of coordinates 
{ X /9 X -, X*} = {20°, 20°, 20°} is used. A singularity 
now exists when a = 0° and p = 0°. 

Clearly, using the X k - 45° port is the preferable 
choice for the off-lateral meridian. If one assumes that 
true angle-of-sideslip limits to be experienced are less 
than ±20°, then for a < 17°, and a > 20°, aroot-closest- 
to-zero criterion is valid for selecting the correct 
root. For values of a between 17° and 20°, the 
Xfy^j = {20°, 20°, 45°} system is indeterminate, and 
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instead the X {jjk] = {20°, 20°, 20°} triple is used. 
With this substitution then, a root-closest-to-zero 
criterion always selects the correct root. 


The root locations for one special case asymmetric 
port arrangement can be analyzed, that is the 
arrangement where Xj =0. Consider the case where, 


<(>. = = cosacosA- ( . 

b i = ±sinX i 

§j = 0, Xj = 0 => dj = cosa 

b j = ° 

- 0 =* a k = cos a cos X k + sin a sin X k 
- cos(a - X k ) 

h = o 

Following the same procedures as before, 


(A-33) 


C 2 
— = cos a 
A 




cos (a-X k ) ^ 2^ 


O',* + r kj ) + r ; , — - rysin x f ) 

cos a 


r*/ sin \) 


2 

cos a 


cos (a -X.) 2 * 

+ r • H—r^sm x,) 

cos a 


r t; -sin 2 X, 


2 

cos a 


cos (a - X k ) 


V cos a 


- 1 


sin 2 A 


(A-36) 


The pressure dependent term, — , can be removed 

1 kj . 

from equation A-34 by substituting in from equation 


C _ 

A] ~ 

_ 2 ~ 2 , i x r- 2 2. (A-34) 

r f .*cos a + T^cos (a- A.*) + r* ; cos a cos A, 

r *; sin % 


Factoring out cos a 


C 2 
— = cos a i 
A 


cos“(a-A,.) 

r ,* + r ; ; 2 + r V 0S 

cos a 


r'kjSin 2 ^ 


(A- 35) 


= cos a 


cos (a - X.) 7 

r « + r ji j— L + r jy (l-$in \) 

cos a 


r*ySin 2 X ( 


Collecting terms, equation A-35 reduces to 


A-5. Defining 


x ¥ = - 


"cos (a - X k ) 

2 ^ 
cos "a 

sin ^7 


-- 1 


then equation A-37 reduces to 


C 

A' 


2 

cos “a 


r r . 

-ilvp-l 

L r */ 


(A-37) 


(A-38) 


and finally 



(A- 39) 


For the X-33 design, there are four possible 
asymmetric off-meridian angle of sideslip triples 
combinations. 
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{[*-;, <D jUK ♦*]} 

= {[20°, ±90°], [0°, 0°], [45°, 0°]} 


(A-40) 


and 


{iKtMi.j.+jUWJ} (A . 41) 

= {[20°, +90°],[0 o , 0°], [20°, 0°]} 

For the first sets of coordinates { [ <|> , 
IKj^jl [&*,♦*]} = {[20°, ±90°], [0°,0°], [45°, 0°]} 
the solutions for P 2 are plotted as a function of p t and 
the local angle of attack in figure A-3(a). In the range of 
ipi<25°, it is interesting to note that the second 
solution is nearly independent of the value of the first 
solution. For 4^ = 90°, P 2 - ~80°, and <J>| = -90°, P 2 = 
80°. 


Figure A-3(b) presents similar results for the p 2 
solutions for the second set of coordinates {[^, <t>J , 
[XjAjl = {[20°,±90 0 U0 0 ,0 0 H20°,0 0 ]}. 

Again, in the range of ip 1 < 25°, it is very interesting to 
note that P 2 is identical to the value produced by the 
previous set of coordinates ( X k = 45°). Clearly, the terms 
involving X k are canceled out of equation A-38. As a 
result, for the asymmetric configuration a root-dosest- 
to-zero criterion will always select the correct solution. 

Angle-of-Sideslip Solution Summary 

The quadratic equation, A-22, has two roots, with no 
guarantees that they be in general orthogonal. However, 
several conclusions can be reached for special geometry 
arrangements which X-33 exploits in its design. These 
conclusions are: 


1. For triples which lie on the lateral meridian, the 
solutions are nearly orthogonal, and for ip I < 25°, 
using a root-closest-to-zero selection criterion will 
always select the correct solution. 

2. For off-meridian triples arrangements, only the 
symmetric arrangement presents a problem. For 
this arrangement, there exists a singularity point at 
which all three incidence angles are equal. In the 
vicinity of this singularity, the system is highly ill- 
conditioned. This ill-conditioning makes the 
system indeterminate and selecting the correct root 
is impossible without prior information. 


3. For the X-33 configuration, using the lower 
45° port to complete the symmetric triple puts the 
singularity at 18.207° local angle of attack. Using 
the lower 20° port puts the singularity at 0° local 
angle of attack. 

4. If one assumes that true angle-of-sideslip limits, to 
be experienced, are less than ±25°, then for the 
X { ^}={20°, 20°, 45°} triple, for a <17°, 
a > 20°, a root-closest-to-zero criterion is valid for 
selecting the correct root. For values of a between 
17° and 20°, the \ ijk} = (20°, 20°, 45°) system 
should be considered indeterminate, and the 

{20°, 20°, 20°} triple should be used 
instead. If this substitution is made, then a root- 
closest-to-zero criterion will always be valid for 
selecting the correct root. 

5. For the asymmetric off-meridian configuration, the 
alternate P solution will always hover around 
±80°. Thus the root-closest-to-zero criterion is 
always valid. 

Static Pressure. Impact Pressure, and Mach Number 
Solution 

Once the local values of a and p have been solved 
for, then the incidence angles at all of the ports can be 
evaluated, and only e, and q c remain as unknowns 
in the pressure equations. Assuming that n pressure 
measurements are to be used, the original model 
(eq. A-l) can be written in matrix form with static and 
impact pressure broken out as variables: 


P\ 


Pn 


^COS~(0J )+^ 
1 

^Esin~(0,) J 



(A -4 2) 


^COS 1 2 (0 ;I )+ > 

v esin 2 (0„) y 


Ignoring for now, the fact that £ is implicitly a 
function of static and impact pressure, equation A-42 can 
be solved for the static and impact pressure using 
weighted linear regression. 19 The result is 
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I <7, "I <?,(«,) 

i = 1 i=l 

-X9.W i W 

i = I ( = I 


X 

L 1=1 


(A-43) 


X* 

Li= 1 J 


X 9/( Q i)' 

Li = ] 


n i2 

i 

L i = 1 


potential instability regions be identified. In general it is 
impossible to analyze the stability of nonlinear equations 
with two unknowns; however, a linearized stability 
analysis will determine the behavior of the system with 
respect to small disturbances. Defining the terms 



(A-46) 


where = cos 2 (0 ? ) + esin (Gy). In equation A-43, 
the q-s are weights which have a nominal value of 1 .0. 
Setting the value of q { to zero, weights the i th pressure 
reading out of the algorithm. Equation A-43 looks like a 
closed form solution, but the £ terms on the right-hand 
side of equation A-43 are implicitly a function of P ^ 
and q c . Thus equation A-43 is actually nonlinear and the 
solutions for P ^ and q c must be extracted iteratively. 
Defining 


and recalling that £ is a function of Mach number, and 
implicitly a function of static and impact pressure, then 


Q, 1 


M[x) = 


(A-47) 


Q 


n 


I 



(Q ( , y) ) 1 


P\ ■ 

. . 0 

M U) = 


, Q = 




(Q^) 1; 


0 . 

• ■ <ln 


where = cos 2 (0 | ) + e ( ^sin 2 (0 f ) , equation A-43 
is written as an iterative estimator 


P. 


Hj + i ) 


M {j)T QM U) 


-i 


M 


U)T, 


p 1 


(A-45) 


Pn 


and equation A-42 can be re-written as 

Z = Af[jr]x (A-48) 

Using the definitions of equation A-46 and A-47, the 
iterative estimation equation can be written in the form 

Z = M[x j) ]x j + n (A-49) 

In equation A-49 ,v (y + /) is the estimate after the 
j + /'th iteration, and M[x {J) ] is the matrix of 
equation A-47 evaluated using the result from the f th 
iteration. Subtracting equation A-49 from A-48, and 
expanding M[x] in a Taylor series about x ^ and 
neglecting terms higher than first order in the 
perturbations, equation A-49 reduces to 


In equation A-45 the subscript (j) refers the result of 
the /th iteration. The matrix A# is defined in 
equation A-44 with £ being evaluated using the values 
for Mach number resulting from the previous iteration 
for P „ and q c . 

Algorithm Stability Analysis 


[VM ( J\ x-x U) )] x + M[x U) ](x-x U + n ) = 0 (A-50) 


Defining the error vectors, 




(A-51) 


Since equation A-45 is to be implemented as part of a 
real-time airdata estimation algorithm, it is essential that 


Equation A-51 becomes the linearized error equation 
for the iteration 
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[VM ( px]x U) + M[x U) ]x U + I) = 0 (A-52) 

and the eigenvalues (roots) of the characteristic 
equation 23 

(A-53) 

will determine the linear stability of system. The 
elements of the Jacobian, V M F x , are evaluated as 


! Mi 


subsonic 


1 + 0.2 


= 1.4Mjl + 0.2M 


2.5 




for(M„<l) 


= 1168.45NT 


2M„- 1 


(A-56) 


supersonic 


7M1- 1 


3.5 


for(A/ ec > 1) 


Substituting equation A-55 into A-54, the Jacobian 
matrix reduces to 




X = 



H c sin e i 




(A-54) 




rae^-i 


J_ 

& 
3 M 


KA-57) 


.„=» (5s22Y— 


) 


^ sin 0 


In equation A-54, the parameter 


de 


is the 


sensitivity of the calibration parameter to Mach number 

and is determined by numerically differentiating the 

calibration parameter (e M ) with respect to Mach 
3M 3M_ 


number. The parameters 


and 


are evaluated 


by differentiating the isentropic flow equation (eq. 9(a)) 
for subsonic flow, or the Rayleigh pitot equation 
(eq. 9(b)) for supersonic flow. For y = 1.4, the resulting 
expressions are 

1 


d <h- 


dM 

dP„ 


°{dMj 


(A-55) 


p*(K) 


where the Mach number derivative is 



sin 20 ,, 


- 3e 0> 

_3M oc _ 


I 

dM 


r^r ■ 

{rj 51 


sin 2 0 


pe (;) ] 1 
dM-\ 


Equation A-57 is dependent only on the free stream 
Mach number and the local incidence angles; thus for a 
given port geometry, the stability of the system 
(eigenvalues equation A-53) is dependent only on the 
free stream Mach number and the local flow incidence 
angles. This conclusion is very important, as it allows the 
stability characteristics of the system to be analyzed 
independently of trajectory. 
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Figure 1. A comparison of the VentureStar™ and X-33 vehicles. 
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Figure 4. Exploded view of X-33 FADS pressure plug design. 
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Figure 5. X-33 FADS computational architecture. 
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Figure 6. Stability of the FADS Mach number iteration. 
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Figure 7. The FADS startup algorithm. 


24 

American Institute of Aeronautics and Astronautics 













(a) Variation with angle of attack. 



(b) Variation with angle of sideslip. 

Figure 10. FADS position error coefficient wind tunnel calibration. 
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(c) Variation with Mach number. 
Figure 10. Concluded. 
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(a) Percentage Mach number error. 



(b) Angle of attack error 
Figure 12. X-33 FADS error estimates. 
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Figure A-3. Comparison of p solutions, off-meridian triples, asymmetric arrangement. 
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